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We study measures of decoherence and thermalization of a quantum system S in the presence of a quantum 
environment (bath) E. The whole system is prepared in a canonical thermal state at a finite temperature. Apply¬ 
ing perturbation theory with respect to the system-environment coupling strength, we find that under common 
Hamiltonian symmetries, up to first order in the coupling strength it is sufficient to consider the uncoupled sys¬ 
tem to predict decoherence and thermalization measures of S. This decoupling allows closed form expressions 
for perturbative expansions for the measures of decoherence and thermalization in terms of the free energies of 
S and of E. Numerical results for both coupled and decoupled systems with up to 40 quantum spins validate 
these findings. 
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Decoherence and thermalization are two basic concepts in 
quantum statistical physics 111]. Decoherence renders a quan¬ 
tum system classical due to the loss of phase coherence of the 
components of a system in a quantum superposition via inter¬ 
action with an environment (or bath). Thermalization drives 
the system to a stationary state, the (micro) canonical ensem¬ 
ble via energy exchange with a thermal bath. As the evolu¬ 
tion of a quantum system is governed by the time-dependent 
Schrodinger equation, it is natural to raise the question how 
classicality could emerge from a pure quantum state. This 
question is becoming more important technologically, for ex¬ 
ample in designing quantum computers |2] where decoher¬ 
ence effects are a major impediment in engineering implemen¬ 
tations. 


Various theoretical and numerical studies have been per¬ 
formed, trying to answer this fundamental question, e.g., the 
microcanonical thermalization of an isolated quantum sys¬ 
tem M, canonical thermalization of a system coupled to 
a (much) larger environment SHI t 3], and of two identical 
quantum systems at different temperatures USB. In our 
earlier work |2^ . we found that at infinite temperature the de¬ 
gree of decoherence of a system S scales with the dimension 
of the environment E if the state of the whole system SA-E 
is randomly chosen from the Hilbert space of the whole sys¬ 
tem. We showed that in the thermodynamic limit, the system 
S decoheres thoroughly. 


In this Letter, we investigate for the first time quantitatively 
how classicality emerges from a pure quantum state when 
5 -f £ is of a finite size and at a finite temperature. We as¬ 
sume that the whole system S-\- E is in a canonical thermal 


state, a pure state at finite inverse temperature j3 E-23], and 
investigate measures of the decoherence and the thermaliza¬ 
tion of S. This canonical thermal state could be the result of a 
thermalization process of the whole system S-\-E coupled to 
a large heat bath, which we do not consider any further. The 
state of the system S is described by the reduced density ma¬ 
trix. The degree of decoherence of the system S is measured 
in terms of d, defined below in terms of measures of the off- 
diagonal components of the reduced density matrix. If d = 0, 
then the system 5 is in a state of full decoherence. The dif¬ 
ference between the diagonal elements of the reduced density 
matrix and the canonical or Gibbs distribution is expressed by 
5. Hence, for a system being in its canonical distribution it 
is expected that both d and 5 are zero. We show that un¬ 
der symmetry transformations that leave the Hamiltonians of 
S and E invariant but reverse the sign of the interaction Hamil¬ 
tonian, conditions which are usually satisfied for example in 
quantum spin systems, the first-order term of the perturbation 
expansion of d^ in terms of the interaction between S and E is 
exactly zero. Therefore it is sufficient to study the uncoupled 
system. We demonstrate that the leading term in the expres¬ 
sions for d^ and 5^ is a product of factors of the free energy 
of E and the free energy of S. Hence, these expressions for d^ 
and 5^ allow one to study the influence of the environment on 
the decoherence and thermalization of S. We study the valid¬ 
ity region of these results obtained from perturbation theory 
by performing large-scale simulations for spin-1/2 systems 
with initial states from a canonical thermal state. 


The Hamiltonian of the whole system 5 -F £ can be ex- 
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pressed as 


H = Hs + He + IHse, ( 1 ) 

where Hs and He are the system and environment Hamilto¬ 
nian, respectively, and XH^e describes the interaction between 
the system S and environment E. Here X denotes the global 
system-environment coupling strength. 

The state of the system S is described by the reduced density 
matrix p = Tr^p , where p = I'P) (‘P| is the density matrix of 
the whole system S + E and Tr^ denotes the trace over the de¬ 
grees of freedom of the environment £. The state of the whole 
system can be written as |'T) = j c{i,p) \i,p) , 

where the set of states {|/,p)} denotes a complete set of or¬ 
thonormal states in some chosen basis, and and De are the 
dimensions of the Hilbert spaces of the system S and environ¬ 
ment E, respectively. We assume that 0$ and De are both 
finite and D = DsDe is the dimension of the Hilbert space of 
the whole system S + E. In terms of the expansion coefficients 
c{i,p), the matrix element {i,j) of the reduced density matrix 
reads Pij = j c*ii,p)cij,p). 

We characterize the degree of decoherence of the system S 

by in 


as a “typical” canonical thermal state Il23l] in the sense that if 
one measures observables, their expectation values agree with 
those obtained from the canonical distribution at the inverse 
temperature j3. In the following we present analytical results 
for (7 and 5 when the whole system S + E is in such a state 

The system S is coupled to an environment E with global 
coupling strength X, and therefore we resort to perturbation 
theory with respect to X. Up to first order in the global system- 
environment coupling strength X, we have 11241] 
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where Hq = Hg + He denotes the Hamiltonian of the uncou¬ 
pled system and bath. If we now consider a unitary transfor¬ 
mation that leaves H^ and He invariant but reverses the sign of 
Hse, then ('T(0)|e-/^" |'T(0)) sa Tre-^^o jo = Z^/D, where 
Zo denotes the partition function of the uncoupled system. The 
expression Eq. © for the canonical thermal state becomes 


mp)) 


V Zo 


( 6 ) 


(7 = 




O5—1 Ds 

E E 


!=1 j=i+l 



Using Eq. I© and Eq. (| 6 l), we find that the first-order term of 
the perturbation expansion in X of the expectation value of ( 7 ^ 
is given by 


where pij is the matrix element {i,j) of the reduced density 
matrix p in the representation that diagonalizes Hs- Clearly, a 
is a global measure for the size of the off-diagonal terms of p. 
If (7 = 0 the system is in a state of full decoherence (relative to 
the representation that diagonalizes Hg). We define a quantity 
measuring the difference between the diagonal elements of p 
and the canonical distribution as dll 


5 = 


\ 



(3) 


where {Ef} denote the eigenvalues of Hs and b is a fitting 
inverse temperature. The quantities C 7 and 5 are respectively 
general measures for the decoherence and thermalization of S. 

A canonical thermal state is a pure state at a finite in¬ 
verse temperature p defined by (the imaginary-time projec¬ 
tion) iiilQ 


|'E(j3)) = / (‘T(0)|e^^"|'E(0))*/^ , (4) 


where |'E(0)) = Y,f=idi\i) and {di} are complex Gaussian 
random numbers. The justification of this definition can 
be seen from the fact that for any quantum observable A 
of the whole system |2j., one has ('T(j3)| A |'T(j3)) 


TrAe^^^/Tre^^^. The error in the approximation is of the 
order of the inverse square root of the Hilbert space size of 
the whole system Kill and therefore the approximation im¬ 
proves for increasing D. One may consider the state |'E(j3)) 


(A'-}), = -!>>■ (g) 


D 


D+1 


(7) 


where Z 5 = Here and in the following ((•)) denotes 

the expectation value with respect to the probability distribu¬ 
tion of the random numbers {dj}. Applying the same unitary 
transformation as discussed before results in ^^( 7 ^))^ = 0 . 
Hence, to study the decoherence of a system S coupled to an 
environment £ up to first order in X it is sufficient to study the 
uncoupled system (A = 0). Eurthermore, one state of the type 
Eq. © is enough to get the correct expectation value provided 
the dimension of the Hilbert space is large iH. 

Therefore, we now focus on the uncoupled system. There 
exist simple relations for the eigenvalues Ej (eigenstates \Ej)) 
of the Hamiltonian H in terms of the eigenvalues Ef, Ep 
(eigenstates |£f), \Ep)) of the system Hamiltonian Hs and 
environment Hamiltonian He, respectively, i.e., Ej = Ef -\-Ep 
and \Ej') = |£?) |Ep ). We first perform a Taylor series expan¬ 
sion of < 7 ^ up to second order in \d'f' about the point 1/D and 
then calculate the expectation value of (7^. A lengthy calcula¬ 
tion gives 

((( 72 ^) = }_e-mFEm-FEm _e-mFs( 2 fi)-Fsm\^ 
p-wiFEm-FEifi)) 

X (e-2m{2P)-Fsm_^-3piFsm-Fs(P))) 
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I -Afi(FE(2fi]-FEm)„-2fi(Fs(2fi]-Fsm 

2{D+\f 

(8) 

where Fs{nfi ) and («j3) represent the free energy of the sys¬ 
tem S and the environment E at the inverse temperature njS, 
respectively. Obviously, in most of the cases the first term 
dominates. 

Two special cases are of interest. First, if j3 = 0, we re¬ 
cover the previous result E{a^) = {Ds —1)/2(D+1) . In 

the vicinity of j3 = 0, the first-order term of the Taylor expan¬ 
sion of Eq. ® vanishes. Hence in the high temperature limit, 
(((j2)) = (Ds-l)/2(D+l) + ^(j32). 

If the temperature approaches zero, Eq. (l8]l becomes 




gs- 1 A_ DsDe 

2-8sgE V [DsDe + 1) gsgE 


(9) 


where gs and gE refer to the degeneracy of the ground state 
of the system S and environment E, respectively. This ex¬ 
pression yields zero if the ground state of the system is non¬ 
degenerate. Eor a system with a highly degenerate ground 
state {gs 3> 1) the expression goes to l/2g£. Eor a system 
with known gj > 1 and a large environment De ^ 1, at small 
X and at low temperatures, if one measures one can 

determine the degeneracy gE of the ground state of the envi¬ 
ronment. This is a new, strong prediction. The ground state 
degeneracy gE of the environment can be esitmated by only 
measuring quantities in S. Eurthermore, Eq. (|9]l should be a 
guide to engineering quantum systems that have limited de¬ 
coherence at low temperatures, namely systems with lar ge g s 
(one characteristic of topological quantum computation 11251] ') 
but non-degenerateg£ so ~ \j2. 

Similarly, we can make the Taylor expansion for 5^ up to 
second order with respect to both \d\^ and b about the points 
1/D and j3. The expectation value of 5^ is given by 

^ ^-2fi{FE{2fi)-EE{fi)) ^^-2fi(Es{2fi)-Es{P)) 

_2e-3j3(F5(3)3)-F5(j3))+e-4j3(f5(2i3)-D(i3))^ 

+e-2PiFs{2li)-Fsm ( 2 j3)/(4kBj3 2 ) 

+ (D5(2j3)-Ds(j3))2)](A/7)2, (10) 

where Ab = b — j5, is the Boltzmann factor, Cs{nj5) and 
Us{nfi) are, respectively, the specific heat and average energy 
of the system S at inverse temperature «j3. It is obvious that 
for the uncoupled system b = j5. Eor the coupled system, b is 
not necessarily equal to j3. 

In order to test the perturbation theory predictions, we nu¬ 
merically simulate a spin-1/2 system divided into a system 
S and an environment £ in a canonical state. We consider 
a quantum spin-1 /2 model defined by the Hamiltonian of 
Eq. ([TJ where 


Hs 


Ns-l Ns 

-LIE WSJ. 


(=1 j=i+la=x.y,z 


( 11 ) 


Ne-1 Ne 

he=-j: e e (12) 

z=l j=i-\-ia=x,y,z 
Ns Ne 

HsE = -tL L N^jSNf- (13) 

1=1 ;=1 cc=x,y,z 


Here, Sf and /“ denote the spin-1/2 operators of the spins 
at site i of the system S and the environment E, respectively. 
The number of spins in S and E are denoted by Ns and Ne, 
respectively. The total number of spins in the whole system 
isN = Ns +Ne- The parameters 7// and D.fj denote the spin- 
spin interactions of the system S and environment E, respec¬ 
tively and AJT denotes the local coupling interactions between 
spins of S and spins of E. 

Erom Eqs. (II 11113b it is clear that the Hamiltonian of the 
whole system, Eq. dD, obeys the symmetry properties which 
are required to make the first-order term of the perturbation 
expansion of the expectation value of (see Eq. (|7]i) ex¬ 
actly zero. Namely, by only reversing the spin components 
of the system or environment spins, Hs and He do not change 
but the sign of Hse reverses. Note that such a symmetry is 
also obeyed in the case that there is no interaction between 
the environment spins, e.g. for an environment Hamiltonian 
He = — I IS 27|. In this particular case, it 

is only required that Hs is an even function and Hse an odd 
function under reversal of all spin components of the system 
spins. 

In our simulations we use the spin-up - spin-down basis 
and use units such that Ti = \ and ke = 1 (hence, all quan¬ 
tities are dimensionless). Numerically, the propagation by 
exp(—)3//) is carried out by means of exact diag onalization 


and the Chebyshev polynomial algorithm Il28l - l32l] with initial 
state |‘F(0)) (see Eq. (|4|i). These algorithms yield results that 
are very accurate (close to machine precision). 


1 /2 

Eigure [T] shows the simulation results for ((( 72 ^^ ' ^nd 

(( 52 /^ obtained by exact diagonalization for the whole 
system S -\- E being a spin chain with Nj = 4 and Ne = 8. 
The system S and environment E consist of two ferromag¬ 
netic spin chains with isotropic spin-spin interaction strengths 
Jf-- = 7 = Elf - = = 1. They are connected by one of their 

end-spins, with an interaction strength A^^ j = A. The global 
system-environment coupling strength X = \. The simulation 
results are averaged over 1000 runs with different initial ran¬ 
dom states. Substituting the numerically obtained values for 
the free energy of the system and environment for XA = 0 in 
the analytical expressions for ((( 72 /y /2 gjygjj 

by Eq. ([^ and Eq. (fl^ . respectively, results in the solid lines 
depicted in Eig.[T] It is clear that the simulation results for the 
uncoupled system (AA = 0) and for the coupled cases when 
(j37) * > 6AA agree with the analytical results for the whole 
range of temperatures. As the temperature decreases, the 
state of the whole system S + E approaches the ground state, 
((( 72 )) becomes constant, its numerical value being given by 
Eq. (01. Eor the case at hand, gs = 5, gE = 9, Ds = 16 and 
De = 256 and Eq. (01 yields ((( 72 ^y /2 „ o.21, in excellent 
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FIG. 1. Simulation results for and ((5^))'^^ for spin- 

1 /2 chains with Ns = 4, Ne = 8, 7 = fl = 1 and various interaction 
strengths A as a function of the temperature l/j3y. The solid lines 
are obtained from Eq. ® and Eq. dlOb by using numerical values for 
the free energies Fs{nj5) and The dotted lines are guides to 

the eye. 


agreement with the numerical data. In the coupled case and 
for small l/jST, ((<7^))*^^ develops a plateau different from 
that of the uncoupled case. The dependence of this plateau on 
j3 or XA is nontrivial, requiring a detailed analysis of how the 
ground state of S + E leads to the reduced density matrix of S 
(in the basis that diagonalizes Hs). In this respect, the p or XA 
dependence of the data show in Fig. [T] are somewhat special 
because the ferromagnetic ground state of the system does not 
depend on XA. 

For the spin system under study with XA^Q, the first-order 
term of the perturbation expansion of the expectation value 
of (7^ in terms of PXA is exactly zero. Hence, for a weakly 
coupled system (AA small) deviations from the analytical re¬ 
sults Eq. ([8]) and Eq. (fTol i obtained for the uncoupled system 
(AA = 0), are, as expected, seen only in the low temperature 
region. The numerical results (symbols) in Eig. 1 are in excel¬ 
lent agreement with the predicted results (solid lines) as long 
as PXA is small. Eor a finite j3A A, the plateaus at low tempera¬ 
ture may or may not be reached and therefore the perturbation 
results are no longer applicable. 

In order to study the behavior of O’ and 5 as a func¬ 
tion of the global coupling interaction strength A, we per¬ 
formed large-scale simulations for a spin system configured 
as a ring with Ns =4- and Ne = 26,36 at the inverse tempera¬ 
ture j317| =6. This required working with a Hilbert space of 
size up to D^Ll X10*^. We assumed that the spin-spin inter¬ 



X 

FIG. 2. Simulation results for o for a ring with Ns = 4, Ne = 26 
(open circles) and Ns = 4, Ne = 36 (solid circles) as a function of 
the global interaction strength A for j5\J\ = 6. For the values of the 
interaction parameters, see text. The solid lines are fits to the data 
as described in the text. The top (bottom) dashed line represents the 
value obtained by simulating the non-interaction system (A = 0) with 
30 (40) spins. 


action strengths of the system S are isotropic, = J = — I, 
that those of the environment E are only nonzero for near¬ 
est neighbors and that only A“^ j and A“^^ are nonzero. The 
nonzero values of njT and AJT are generated uniformly at ran¬ 
dom from the range [—4/3,4/3]. 


In Fig. ID we present the simulation results for a as a func¬ 
tion of A. Least square fitting of the data for <7^ to polyno¬ 
mials in A, we find that a polynomial of degree 7th y ields the 


best fit, for both the 30- and 40-spin system data ||33LI34|] . The 


behavior of 8 is very similar to that of a and is therefore not 
shown. From Fig. 2 it is clear that for A Ri 1, the scaling of 
(7 with the dimension of the Hilbert space of the environment 
is almost completely suppressed. This is a finite temperature 
effect: for j3 = 0, this scaling property holds independent of 
the coupling A 


Summarizing, in this Letter, we investigate the measures 
of decoherence and thermalization of a quantum system S 
coupled to an environment E at finite temperature. If the 
whole system S + E is prepared in a canonical thermal state, 
we show by means of perturbation theory that <7^, the de¬ 
gree of the decoherence of S, is of the order jS^A^. Up to 
the first order in the system-environment interaction we find 
(7^, 5^ exp{—2j3 [Fe(2P) — Fe{P)]} where Fe is the free en¬ 
ergy of the environment. This provides a measure for how 
well a weakly-coupled specific finite environment can deco¬ 
here and thermalize a system at an inverse temperature j3. A 
measure for how difficult it is to decohere a quantum system 
is given by ratios of free energies of the system, and at low 
temperatures by <7^ ^ 1 /2gE for a highly degenerate ground 
state for S and a ground state degeneracy gE for E. We per¬ 
formed numerical simulations of spin-1/2 systems to assess 
the region of validity of the analytical results. 


Strictly speaking, the system S completely decoheres if 
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there is no interaction between S and E and if Ne —t- If 
S is coupled to E, the S-E interaction is important. Generally, 
<7 and 5 for a finite system S are finite even in the thermody¬ 
namic limit (Ne °°)- However, if the canonical ensemble 
is a good approximation for the state of the system for some 
inverse temperatures j3 up to some chosen maximum energy 
Ehoid > 0 (measured from the ground state), then it is required 
that exp(— jSii/joM) ^ <7. By determining the crossover of the 
left- and right-side functions, we find a threshold for the tem¬ 
perature above which the state of the system is well approx¬ 
imated by a canonical ensemble, and below which quantum 
coherence of the system is well preserved. 
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